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I. INTRODUCTION steepness, and the first criteria of Michell, J.H. (1893). 
Wavelength?) sad waveadmberk inthe ke an The results of ie research were widely eed by 
L other researchers in developing breaking criteria in 

i " shallow waters, including Miche (1944), Battjes and 
accuracy of water wave analysis is largely determined by Jansen (1978), Ostendorf and Madsen (1979), Battjes and 


these parameters. The equation for calculating these Stive (1985) and Rattanapittikon and Shibayama (2000). ). 
parameters is called the dispersion equation and the most 


: A : : The newest criteria of critical wave steepness is the criteria 
used is the dispersion equation from the linear wave theory f . 
(Dean (1991)) proposed by Toffoli, A., Babanin, A., Onorato and 


Waseda,T. (2010), with the same form as Michell's (1893) 
criteria, but having different coefficients. In this study, the 
criteria from Toffoli et all will be used. 


a very important water wave parameter, where the 


Hutahaean (2019) formulated the dispersion equation for 
water waves using the weighted total acceleration in the 
KFSBC equation and the Euler momentum equation, 
which the weighted total acceleration in both equations has 
the same form. In this study, the weighted total 
acceleration was developed in different forms. 

With the KFSBC equation and Euler's momentum 
equation using weighted total acceleration, the dispersion 
equation for water waves was formulated. To get the 
weighted coefficient value, the critical wave steepness of 
the equation is calibrated against that of the previous 
research results. There are two critical criteria for wave 


Wiegel (1949,1964) proposed an equation for the 
maximum wave height for a wave period. By using this 
wave height maximum, wavelengths that meet the criteria 
of Toffoli, et al (2010) can be calculated. So that we get a 
wavelength value that meets the maximum wave height 
criteria and the critical wave steepness criteria. 
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II. MAXIMUM WAVE HEIGHT AND CRITICAL 
WAVE STEEPNESS 

a. Wiegel’s criteria (1949,1964) 

According to Weigel (1949,1964) , the maximum wave 


height in a wave period is, 
gT? 
15.62 
T = wave period (sec) 


Hmax = 
g = gravitational force (m/sec?) 


b. Michell’s criteria (1893) 
According to Michell (1893) critical wave steepness is 


c. Toffoli, A., Babanin, A., Onorato and Waseda,T. 
(2010) 

According to Toffoli, A., Babanin, A., Onorato and 
Waseda,T. (2010), critical wave steepness is 


H 0.170 ....... (3) 


-= 
L= wavelength 

With the input wave height from (1), it can be calculated 
the wavelength with (2) or (3), it is then compared with the 
deep water wavelength of the dispersion equation from the 
linear wave theory (Dean (1991)), 


T2 
Ly = o—...(4) 
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(1). It can be seen that the wavelength produced by the 
linear wave theory is too long, so that the critical wave 
steepness criteria of both Michel and Toffoli will never be 
achieved. 


Furthermore, if the wavelength of (4) is used to calculate 
the critical wave height with (2) and (3), a very large 
critical wave height is obtained as presented in table (2) 
below. 


Table.2: Wave height maximum according to (2) and (3), 
wave length is calculated by (4) 









































T Hmuicn Hrof 
(sec) (m) (m) 
6 7.98 9.56 
7 10.86 13.01 
8 14.19 16.99 
9 17.96 21.5 
10 22.17 26.54 
11 26.83 32.12 
12 31.93 38.22 
13 37.47 44.86 
14 43.45 52.02 
15 49.88 59.72 











Table.1: The wavelength from (2), (3) and (4) 















































T Hmax | Lrof | Lmicn Lo 
(sec) (m) (m) (m) (m) 
6 1.45 8.54 10.22 56.21 
7 1.98 11.62 13.91 76.5 
8 2.58 15.18 18.17 99.92 
9 3.27 19.21 22.99 126.47 
10 4.03 23.71 28.39 156.13 
11 4.88 28.69 34.35 188.92 
12 5.8 34.15 40.88 | 224.83 
13 6.81 40.07 47.98 | 263.86 
14 7.9 46.48 55.64 | 306.02 
15 9.07 53.35 63.87 | 351.29 











On table (1), Lror is the wavelength that is calculated by 
(3), Lmicn Wave length that is calculated by using (2), and 
Lo is the wavelength that is calculated by using (4). The 
results of the wavelength calculation are presented in table 
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In table (2), it can be seen that the maximum wave height 
that is calculated using the wavelength of the linear wave 
theory produces a very large wave height. In table (2), 
Hmich the wave height is calculated using the Michell’s 
criteria, while the wave height Hro¢ is calculated using the 
Toffoli’s criteria. Considering that Toffoli et al's criteria 
produce a greater maximum wave height than Michell's 
criteria, this research uses the Toffoli et al’s critical wave 
steepness criteria to calibrate the wavelengths produced by 
the dispersion equations of this study. 


Il. WEIGHTED TOTAL ACCELERATION. 
Courant (1928) stated that the relation between 6x 
and ôt in the water wave equation is ôx = 3.0 C dt where 
C is wave celerity, x is the horizontal axis and t is the 
time. 


Hutahaean (2021) who used the Taylor series on a 
sinusoidal function 

f(t,x) = A cos ot cos kx ....... (5) 

A= amplitude 

o = angular frequency= =. T = period 


k = wave number 
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gained a relationship that is more or less the same as the 
Courrant relationship, that is 6x = 3.01519 C ôt 


Hutahaean (2021) used the Taylor series for functions that 
take the form of: 

f(t,x,z) = A cos ot cos kx cosh k(h + z) ...... (6) 

which is the solution to Laplace's equation (Dean 1991), 
where z is the vertical axis, get the relation or 6z = 
3.0012 C dx or 6z = 3.012? C ôt. 

Based on Courrant (1928) and Hutahaean (2021), the 
relationship obtained 

ôx = y C ot... (7) 

ôz = y? C ôt....(8) 


where y is called as weighting coefficient. 


Relationship between ôx and dt and between ôz and ôt 
were done on Taylor series order 1, in a function of time 
and space f (t, x) by working on the weighting coefficient 
on time differential terms. 


(t +ôt,x + 6x) = f(t,x) + TERT. af 
f a a er On ae 
While on a function f(t, x, z) 


f(t + 6t,x+ 6x,z+ 6z) = f(t,x,z) + 





af af df 
2e 2 Eea Z 
y ôt r + yx + ÔZ a7 
a. Function F = f (t,x) 
df df 
F(t + St, x + dx) = f (t,x) + yot Te + dx 
f(t+dt,x+6x)-f(t,x) af xaf 
bt ~ ae" Stax 
i ae an (9) 
dt dt dx 


For the water level equation n = n(x, t) 
Dn dy % dn 
dt Vat “ax 





KFSBC, Dean (1991) is : 

w= a+ Uy = meer (10) 

Where w,, is the the velocity of water is in the direction of 
the vertical axis-z - on the surface of the water, whereas Un 
is the velocity of water particles in the direction of the 
horizontal axis-x on the surface of the water. It can be seen 
that KFSBC is the total derivative in the Taylor series for a 
functionF = f(t,x) as in (9). Therefore, by using the 
weighting coefficient as in (9), KFSBC becomes 


d d 
Wy = V yet Un se mee (11) 


b. Function F = f (t,x, z) 
Taylor series firstorder for function f (t, x, Z) 
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f(t + 6t,x+ 6x,z+ 6z) = f(t,x,z) 
20, df df df 
+y* ot J + yOx + ôzZ o re (12) 





The total acceleration is, 

bf _ af, af, af 

dt Y at ax az 

The total acceleration of the water particles in the 


horizontal direction-x is, 


In the same way, the total weighted total acceleration can 
be obtained in the vertical direction.-z 

Dw dw dw dw 

as y? atg Wel) 

With (13) and (14), then the equation for Euler's 
momentum in the horizontal direction-x and the vertical 
direction-z becomes, 








2 du du du _ _1dp 

y at eg tT Wa oie ooe (15) 
2 dw dw dw _ idp | 

y atutu" E O a (16) 


Equations (15) and (16) are a modification of Euler's 
momentum equation in the total acceleration term, by 
adding the weighted coefficient y. The difference between 
this study and Hutahaean’s (2019a) is that in this study the 
weighting coefficient is not only in time differential terms. 


IV. FORMULATION OF SURFACE MOMENTUM 


EQUATIONS 
The irrotational flow condition is substituted in (15) 
and (16) where, E ; 
dz dx 

amuy hed =, b82 
a Ton (yuu + ww) = a (17) 

ziw gt ded ze 
Y +55, (yuu + ww) = oa I (18) 


Eq. (18) is written into an equation for p and integrated 
with respect to the vertical axis z from the elevation z to 
the water surface 7, and the dynamic boundary condition 
of the surface where the surface pressure p, =0 is 


applied, thenthe equation of pressure p is obtained 


p "dw 1 
5 =y? l a + z unun + Wy Wy) 


1 
—5 (yuu + ww) +g- z) 


This equation is differentiated with respect to the 
horizontal axis-x, 
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1dp id [Mar 
Z 


pax Y dx dt 


1d 
+ 24x (yunun + WyWp ) 


-22 (yuu +ww)+ g% on (19) 
Substituting (19) by (17), 
,du_ sd ("dw 
Pi ae) ae 
1d 


a 

-zg (YUnun + Wy Wp) = g (20) 

The solution to f 7 SY dz was done using the equation of 
dx°Z dt 

the potential velocity of water waves (Dean (1991)), that is 

@(x,z,t) = Gcoshk(h+z)coskxsinot _....... (21) 

The velocity and acceleration in the vertical direction-z 

are, 


w=- = = —Gksinhk(h + z)coskxsinot..(22) 


dw 


a” —Gksinhk(h + z)øcoskxcosot ........ (23) 


—G(coshk (h +n) — coshk(h + z))coskxcosot 

then was differentiated with respect to the horizontal axis- 
x, 

d ("dw 
dx J, dt 


Gk(coshk(h +n) — coshk(h + z)) sinkx cosot 


Horizontal direction velocity and acceleration-x is, 

u= -— z = Gkcoshk(h + z)sinkxsinot.....(25) 

= Gkocoshk(h + z)sinkxcosot.....(26) 

From(24) and (26) the following relationship is obtained 


d d du d 
pP Zaz =a sn (27) 
dx~Z dt dt dt 


Substituting (27) to (20) the equation of 
momentum is obtained, 


surface 


2 dun 1d 


= dy 
a + agg (Uy + wawy) = -g Z-(28) 


V. FORMULATION OF THE DISPERSION 
EQUATION 


The dispersion equation is formulated at the 





characteristic point, where cos kx = sin kx = cos ot = 
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; v2 : ae ; 
sinot =—. At this characteristic point, water surface 
f A 
elevation 7 = P 


The dispersion equation will be formulated using two 
equations, namely the KFSBC and the surface momentum 
equation, (28). 

a. KFSBC operation. 


Hutahaean (2021), by integrating KFSBC (11) against time 
t, obtained 


n(x, t) = Acoskxcosot...(29) 


At the characteristic point, 
an _ = kA 

yes .(30) 

Where, 


A = = coshk (n+<)(1- 5 oor G1) 
Substituting (31) to (30) 
an — ESE coshk (h + 2) (1 = “*) = (32) 


dx 2 oy 


Equation (31) can be written as an equation for G, 
20yA 


G= r EO) 


~ kcoshk(h+5)(1- * 


b. Operation of Surface Momentum Equations 


Substituting (25) for horizontal velocity u, (22) for vertical 


velocityw and (32) for to (28), where at the 


. e . A 
characteristic point 7 = > 





ip (ier nk (n+) +5 ae 
yrotr=Y- cos =] + -—— 
2 2 2 coshk (h + 2) 
_ gk (: J] 
~ 2oy 2 
Substituting G by (33), 
kA 
y2o(1- —) + ovAk = 1 + — oa 
2 cosh*k (h + 2) 
_ gk (: ), 
~ 20y 2 
ak(n+4) 
Hutahaean (2021) found that 2°=0, so the 


Ox 


A : 
coshk (h + 2) values is constant. In deep water, element 


1 : ; ; 
= 0, is also constant. The final equation 


cosh?k(h+5) 
becomes, 
2 (1 a+ noe (1 a 
yo 2) OIYAY A E Soy 2 
Pee (34) 
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This equation is the dispersion equation for deep water, 
which can be solved by the Newton-Rhapson method. The 
solution using the Newton-Rhapson method requires an 
initial estimate value of k. The initial estimation value can 
be obtained by working with the assumption that in deep 
water the convective acceleration term, the second term on 
the left side (34) can be ignored, so that a simpler equation 
is obtained, 

2y30* = gk (1 — =) 


Or, 
oe — gk+2y3o? = 0 ......(35) 
Calculating k using (35), when the formulas a, b, c is used 
are: 

-b -Vd 
k = —— 

2a 

The determinant d in (35) is: 
d = g? —4gAy3o? 
At the determinant value d = 0, the maximum amplitude 
for a wave period will be obtained, that is 


Amax = mE ae (36) 


VI. ADJUSTMENT TO CRITICAL WAVE 
STEEPNESS CRITERIA AND MAXIMUM 
WAVE HEIGHT IN DEEP WATER. 


In (34) there is a weighting coefficient y whose 
value will determine the resulting wavelength. In this 
section the weighting coefficient value y will be 
determined using the maximum wave height criteria and 
critical wave steepness criteria. 


The calculation was done by trial and error that is by 
determining y, the wave amplitude A is calculated using 
(36), then the wave number kis calculated using (34), 


where the wavelength L = = then the wave height H is 


calculated using (37)-(40). Trial and error was carried out 
until wave height H = Hyg, was obtained, whereHmax 
was calculated by using (1). 


Calculation of wave height H was done using the water 
level equation 


No(x,t) = Agcoskxcosot....(37) 


Ag is the wave amplitude calculated by (36). Then the 
wave amplitude A was calculated, 


Water wave surface elevation is: 
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n(x,t) = Acoskxcosot....(39) 


The calculation was done fort = 0. By using the water 
wave surface equation the 7-maximum and 7-minimum 
were calculated. Wave height H is, 

H = Nmax — Nmin----(40) 

The trial error on value y was done until H = H,,,, where 
the condition was achieved at y = 1.4614, but the critical 


wave steepness == 0.210 is greater than the Toffoli et 


all’s (2010) criteria that is = = 0.170. 


Table.3: Wave height and critical wave steepness ony = 
1.4614 


T L A H | Hmax 





H 


(sec) | (m) m) | (m) | (m) L 
6 6.894 | 0.717 | 1.451 | 1.451 | 0.21 
7 9.383 | 0.975 | 1.975 | 1.975 | 0.21 
8 12.256 | 1.274 | 2.58 | 2.58 | 0.21 
9 15.511 | 1.612 | 3.265 | 3.265 | 0.21 
10 19.149 | 1.99 | 4.031 | 4.031 | 0.21 
11 23.171 | 2.408 | 4.877 | 4.878 | 0.21 
12 | 27.575 | 2.866 | 5.804 | 5.805 | 0.21 
13 | 32.362 | 3.364 | 6.812 | 6.812 | 0.21 
14 | 37.533 | 3.901 7.9 | 7.901 | 0.21 
15 | 43.086 | 4.478 | 9.069 | 9.07 | 0.21 
























































The water wave surface profiles was calculated by (37), 
(38) and (39) are sinusoidal in shape where“ = 2, but with 


truncated peaks. This condition is in accordance with the 
characteristics of the sinusoidal solution to Laplace's 
equation (Dean (1991)). 


15 

1 

0.5 

£ o 
= 

-0.5 

-1 


-1.5 
x (m) 


Fig.l: Wave profile for T = § sec., y = 1.4614 
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In this following section, trial-error of the weighting 
coefficient value was done until wave steepness value 


Pinas — 0.170, which is the Toffoli criteria. This condition 


was achieved atthe weighting coefficient value y = 1.45, 
but wave height H is less than Hmay, as can be seen in 
Table (4). 


Table.4: Wave height and critical wave steepness on y = 
1.45 


T L A H | Hmax 





Amax 


(sec) | (m) m) | (m) | (mM) L 
6 7.101 | 0.734 | 1.206 | 1.451 | 0.17 
7 9.665 | 0.998 | 1.642 | 1.975 | 0.17 
8 12.623 | 1.304 | 2.144 | 2.58 | 0.17 
9 15.976 | 1.651 | 2.714 | 3.265 | 0.17 
10 | 19.724 | 2.038 | 3.35 | 4.031 | 0.17 
11 | 23.866 | 2.466 | 4.054 | 4.878 | 0.17 
12 | 28.403 | 2.934 | 4.825 | 5.805 | 0.17 
13 | 33.334 | 3.444 | 5.662 | 6.812 | 0.17 
14 | 38.659 | 3.994 | 6.567 | 7.901 | 0.17 
15 | 44.379 | 4.585 | 7.538 | 9.07 | 0.17 
























































It can be concluded thatstudy on deep water results 
weighting coefficient valuey in the range 1.45-1.464, with 
a critical wave steepness in the range 0.17-0.21. Where 
with this weighting coefficient value the wave height 
maximum is in accordance with the wave height 
maximum criteria of Wiegel (1949,1964) and the wave 
steepness is in accordance with the Toffoli’s (2010) 
criteria . 


VII. THE EFFECT OF WAVE AMPLITUDE ON 
WAVELENGTH 


In the previous section, it has been found that the 
wavelength and wave amplitude are determined by the 
weighting coefficient value y. In (34), there is an 
interaction between wave ampltudeAand wavelengthL. In 
the following section, the effect of wave amplitude on 
wavelength is investigated. This is important, considering 
that (34) there is no water depth parameter h. So it is 
necessary to explain how the effect of water depth h on 
wavelength L. 
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12 
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0 0.5 1 15 
A (m) 


Fig. 2: Wavelength graphic L towards the wave amplitude 
A, wave perod 8 sec., y = 1.45. 


The effect of wave amplitude on wavelength is shown in 
Fig. (2) where the larger the wave amplitude, the shorter 
the wavelength. When waves move from deeper waters to 
shallower waters, an enlargement of the wave height and 
shortening of the wavelength occurs, known as the 
shoaling phenomenon. Thus the effect of water depth on 
wavelength is through the wave amplitude, where the 
wave amplitude is influenced by the water depth. In other 
words, to obtain wavelengths in shallow water, shoaling 
analysis must be carried out. However, the wavelength 
calculation can be done in a simple way by ignoring the 
wave amplitude enlargement, which will be discussed in 
section 8. 


Table.5: Wave height and critical wavelength on y = 1.00 
T L A H Hmax 





Amax 
(sec) | m) (m) m) | (m) L 


6 28.104 | 2.236 | 4.027 | 1.451 | 0.143 
7 38.252 | 3.044 | 5.481 | 1.975 | 0.143 
8 49.962 | 3.976 | 7.159 | 2.58 | 0.143 
9 63.233 | 5.032 | 9.061 | 3.265 | 0.143 
10 78.066 | 6.212 | 11.187 | 4.031 | 0.143 
11 94.459 | 7.517 | 13.536 | 4.878 | 0.143 
12 | 112.414 | 8.946 | 16.109 | 5.805 | 0.143 
13 | 131.931 | 10.499 | 18.905 | 6.812 | 0.143 
14 | 153.008 | 12.176 | 21.926 | 7.901 | 0.143 
15 | 175.647 | 13.977 | 25.17 | 9.07 | 0.143 



























































To further clarify the effect of wave amplitude on wave 
length, the calculation is done with the weighting 
y= 1.00. With this weighting 
coefficient value, the effect of the convective acceleration 


coefficient value 


Page | 214 


Syawaluddin Hutahaean 


in (34) is neglected with the calculation results presented 
in Table (5). It can be seen that even though it is used y = 
1.00, the wavelength of (34) is still much shorter than the 
wavelength of the linear wave theory (4), where the 
wavelength of the equation can be seen in Table (1). This 
shows the magnitude of the wave amplitude influence on 
wavelength. With this weighting coefficient value, a very 
large wave height is produced, much bigger than Hmax 
Wiegel's criteria (1949,1964), but with a critical wave 
steepness that is very close to Michell's (1893) criteria, 


namely = = 0.142. 


VIII. SIMPLE METHODS FOR CALCULATING 
WAVELENGTHS ON SHALLOW WATER 
The dispersion equation obtained was the equation 
for deep water. In the previous section, it has been found 
that the effect of waterdepth on wavelength is through 
wave amplitude, where the calculation of wavelength in 
shallow waters must go through shoaling analysis. In this 
section analysisof wavelength was done by ignoring the 
phenomenon of wave amplitude enlargement, or by 
working on constant wave amplitudes. The calculation was 
done using the wave number conservation equation 
(Hutahaean (2021)), namely, 
ak(n+2 
au = 0.....(41) 
For waves moving from water depth h, to shallower water 
depth h, this applies 


A, A, 
kp (hz + Z) =k (i, + 3) 


14 
12 


m) 


10 


wavelength L ( 


8 
6 
4 
2 
0 

0 5 10 15 


water depth h (m) 


Fig.3: Wavelength graphic on water depth, T = 8 sec., 
y = 1.45 


By ignoring the change in wave amplitude 
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In (43) it can be seen that forh, less than h, there will be 
an enlargement of the wave number k or shortening on 
wavelength L. 

On Fig. (3), it presents the calculation results with (42) for 
waves with wave period T = 8 sec., weighting coefficient 
y = 1.45, with a linear change in water depth, it is found 
that the wavelength also changes linearly. 


IX. CONCLUSION 
The main conclusion is that the total acceleration 
equation in both the KFSBC equation and the momentum 
equation, which uses the weighting coefficient produces a 
dispersion equation that can produce a wavelength that is 
close enough to the criterion of critical wave steepness. 


The weighting coefficient value obtained was 1.45-1.464, 
where with this weighting coefficient value the critical 
wave steepness and maximum wave height were obtained 
according to the existing criteria of previous researchers. 


In the dispersion equation obtained, there is no effect of 
water depth, but there is an influence of wave amplitude. 
The effect of wave amplitude on wavelength is quite 
significant, that even though a weighting coefficient of 1.0 
is used, the wavelength is still much shorter than the 
wavelength of the linear wave theory. 


The effect of wave amplitude on wavelength is that the 
bigger the wave amplitude is, the shorter the wave length. 
This is consistent with the shoaling-breaking phenomenon 
in waves moving from a deeper water depth to a shallower 
water depth, where an enlargement of the wave amplitude 
occurs which results in a shortening of wavelengths which 
eventually breaking occurs when the critical wave 
steepness is exceeded. 


To get a more definite weighting coefficient value, it is 
necessary to examine the shoaling-breaking phenomenon 
in a model developed using the KFSBC equations and 
Euler's momentum using weighting total acceleration and 
comparing the model results with the breaker index 
equation. 
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